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We investigate four-quark systems consisting of two heavy anti-bottom 
quarks and two light up/down quarks. We propose to solve a coupled 
Schrodinger equation for the anti-bottom-anti-bottom separation using po¬ 
tentials computed via lattice QCD in the limit of static anti-bottom quarks. 
This coupled Schrodinger equation allows to incorporate effects due to the 
heavy anti-bottom spins. First exploratory numerical tests are discussed. 
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1. Introduction 

In recent papers [1, 2, 3, 4, 5] we have studied heavy tetraquark can¬ 
didates combining lattice QCD and quark model techniques proceeding in 
two steps. First we have computed potentials of two static antiquarks QQ 
in the presence of two quarks of finite mass qq (q € {u, d} throughout this 
work) using lattice QCD [1, 2] (such potentials have also been computed by 
other groups, cf. e.g. [6, 7, 8, 9, 10, 11, 12]). The static approximation is 
expected to be a rather good approximation for QQ = bb and allows for a 
comparably easy computation of the potentials. For larger QQ separations 
some of these potentials can be interpreted as potentials of two B and/or 
B* mesons, which are degenerate in the static limit. In a second step we 
have inserted these potentials into the Schrodinger equation for the relative 

* Presented at Excited QCD 2015: Tatranska Lomnica, Slovakia. 
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coordinate of the two B/B* mesons. We have then checked, whether they 
are sufficiently attractive to host bound states, which would indicate stable 
qqbb tetraquarks. For a specific potential (isospin 1 = 0 , light quark spin 
j = 0 ) we have found a bound state with confidence level ~ 2 a and binding 
energy ~ 90MeV, while there seems to be no bound state in any of the 
other channels [3, 5]. 

In the static limit effects due to the spin of b quarks are neglected, e.g. 
there is no mass difference of B and B* for infinitely heavy b quarks. These 
effects, however, could be of the same order as the ~ 90 MeV binding energy 
of the tetraquark predicted in [3, 5], as can e.g. be estimated from the mass 
difference djb* - mg ~ 50 MeV. The goal of this work is to take the heavy 
b spins into account, in particular to estimate their effect on the binding 
energy of the above mentioned (I = 0 , j = 0 ) qqbb tetraquark. 


2. Incorporating heavy bb spin effects 

2.1. Relating qqQQ potential and creation operators 

Due to static quark symmetries it is essential to couple the light spin 
indices and the static spin indices separately, when defining qqQQ potential 
creation operators. To interpret the meson-meson structure generated by 
such operators, one needs to express them in terms of static-light bilinears. 
We do this by using the Fierz identity 


LabScd (Qcin^A (rif) (Qd^Wb fa)) 

= ^Tr (r 6 5 T r^L) (Q^rVVi)) (Q(^)rV 2) (r 2 )), 


( 1 ) 


where A, B, C, D denote spin indices, 

F a € {75,7075,1,7o,7j)7o7j.7j75,7o7j75} and T a is the inverse of T a l . 
The left hand side of this equation has the structure of a qqQQ potential 
creation operator (cf. e.g. [5], eq. (6)), while the right-hand side allows to 
read off, which linear combination of B meson pairs it excites. 

In the following we are interested in those matrices L and S generating 
B and/or B* mesons (Q( 1 + 70)759 and Q( 1 + 70 ) 7 ?+ respectively). After 
some linear algebra one finds that there are 16 such combinations, L, CS T C G 
{C( 1 + 70 ) 75 ,C (1 + 70 ) 7 . 7 } (C denotes the charge conjugation matrix). The 
qqQQ potentials, which have been computed in the static limit, depend only 
on the light spin coupling L, but not on the heavy spin coupling S. There 
are two different potentials, (1) V 5 (r) (corresponding to L = C( 1 + 70 ) 75 ), 
attractive for isospin 1 = 0, repulsive for isospin 1 = 1, and (2) Vj(r) 

1 Similar techniques have recently been applied to relate meson-meson and diquark- 
antidiquark creation operators [13]. 
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(corresponding to L = C(1 + 70 ) 7 j), repulsive for isospin 1 = 0 , attractive 
for isospin 1=1, where r = \ f\ — r 2 1 - 

Note that it is not possible to choose S and L in a way that exclusively 
B mesons appear on the right hand side of eq. (1). One always finds linear 
combinations of B and B* mesons, e.g. for L = CS T C = C( 1 + 70)75 the 
right hand side of eq. ( 1 ) is proportional to R(ri).B(r 2 ) + B*{r\)B*{f 2 ) + 
B*(fi)B*(r 2 ) + B*(fi)B*(r 2 ) (the indices x, y, z denote the spin orientation 
of B*). Taking this mixing of B and B* mesons into account, which differ 
in mass by ~ 50 MeV, is the goal of this work, as already mentioned in the 
introduction. 


2.2. The coupled channel Schrodinger equation 
We study a coupled channel Schrodinger equation 

H^(ri,r 2 ) = E^(ri,r 2 ), ( 2 ) 

where the Hamiltonian H acts on a 16-component wave function \k. The 
components of 'll correspond to the 16 possibilities to combine 
(B (n) , B* (rq), B*(fi),B*(fi)) and (R(r 2 ) , B* (r 2 ),B*(r 2 ), B* z (f 2 )), he. the 
first component corresponds to H(ri)H(r 2 ), the second to H(ri)H*(r 2 ), etc. 

The Hamiltonian can be split in a free and an interacting part, H = 
Hq + H[ nt . The free part is given by 

H 0 = M <g) 1 + 1 < 8 ) M + ^ (m ® l) 1 + y (l (g) m) 1 (3) 

with M = diag(m b, ms*, ms* The interacting part can be written 

according to 

H int = T"V(r)T, (4) 

where 

V(r) = diag (v b {r),...V b (r),V j {r),...V j {r)\ (5) 

V '-V-' v" ' 

4x 12 x 

and T is a 16 x 16 matrix relating the 16 choices for L,S (cf. section 2.1 
and eq. (1)) to the 16 components of T (the entries of T can be computed 
using the Fierz identity (1)). 

3. Numerical solution of the coupled channel Schrodinger 

equation 

Rotational symmetry allows to bring the coupled channel Schrodinger 
equation ( 2 ) to block diagonal form, i.e. to split it into independent simpler 
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equations corresponding to definite total spin J and isospin I 2 


• a single 2 x 2 coupled channel equation: 

J = 0,7=1, meson pairs B(f\ )B{v 2 ) and B*(r\)B* (rq); 

• three identical lxl, i.e. uncoupled equations: 

J = 1, I = 1, meson pairs B(r\)B*(r 2 ) and B*{r\)B{r 2 )] 

• three identical 2 x 2 coupled channel equations: 

J = 1, / = 0, meson pairs 5 (ri)B*(r* 2 ), B*{f\)B{r 2 ) and B*{fi)B*{f 2 )\ 

• five identical lxl, i.e. uncoupled equations: 

J = 2, / = 1, meson pairs B*{f{)B*{f 2 ) ■ 


For the remainder of this section we focus on the J = 0 coupled channel 
equation, where 


Ho,j=o - 

Hint,J=0 


2 m B 0 \ (Pi 

0 2 ms* J y 2 

/(l/4)(y 5 (r) + 3^(r)) 
V(V3/4)(F 5 (r)-VS(r)) 


(1/rriB 0 \ 

2 J V 0 V^b* y ’ 

(\/3/4)(F5(r)-V-(r))\ 

(l/4)(3 Vstf + Vjir))) 


( 6 ) 

( 7 ) 


Introducing center of mass and relative coordinates the partial differential 
equation in fq and F 2 can analytically be reduced to an ordinary differential 
equation for r, 




2m b* 



+ 7/int,J=0 


X{r) = Ex(r ), 


( 8 ) 


where the first component of y represents a B(r\ )B(r 2 ) pair and the second 
component a B*{f\)B*{F 2 ) pair. Following standard textbooks on quantum 
mechanics one can show that the radial wave function of an s wave bound 
state is subject to the boundary conditions 


X(r) ~ 



as r —>• 0 


lim x( r ) 

r —>00 



(9) 


with A, B G R. 

First exploratory numerical tests have been performed with 1 = 0 poten¬ 
tials, i.e. with an attractive V 5 and a weakly repulsive Vj 3 . We integrate 

2 J and I are related, because quarks are fermions and have to obey the Pauli principle 
(cf. [5] for a detailed discussion of quantum numbers of qqbb tetraquarks). 

3 Even though this (J = 0 ,1 = 0) channel is excluded by the Pauli principle, it is 
conceptually interesting to compare numerical results with existing results from [3, 5], 
where heavy spin effects have not been taken into account. 
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eq. ( 8 ) using the Runge-Kutta-Fehlberg method starting with the linear 
asymptotic behavior (9) at tiny r = e>0tor = r max with sufficiently 
large r max >10fm. This integration is iterated many times as part of a 
standard shooting procedure to find parameters A/B and E such that also 
Xl (T max) = Xl(^max) = 0 is fulfilled. 

In Fig. 1 we show results obtained with an unphysically strong attractive 
potential V 5 (roughly a factor 1.5 stronger than the lattice QCD result for 
V 5 ). The intersection of the red line and the green line at E ~ 10.4 GeV cor¬ 
responds to Xi( r max) = Xi( r max) = 0, i.e. represents an energy eigenstate. 
Since E < 2 ms ~ 10.6MeV (2ms is the upper boundary of the plot), this 
eigenstate is a bound four quark state. 
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Fig. 1. Isolines X'i(r max ) = 0 (red) and X2(r ma x) = 0 (green) in the A/B-E plane 
(for unphysically strong attractive potential V 5 ). 

When repeating these calculations with crude fits to the lattice QCD 
results for V 5 and Vj, the situation is less clear, i.e. E ~ 2mg. A more 
careful analysis and treatment of statistical and systematic errors similar to 
what has been done in [5] is needed, to confirm or rule out a bound state. 
In any case, one can conclude that the heavy bb spins counteract four-quark 
binding. 




. . 

.. 


Xl ( r max) ~ ^ " 

X2( r max) ~ ® 




4. Outlook 

Most interesting will, of course, be an investigation of the physical chan¬ 
nels listed at the beginning of section 3, in particular the (J = 1,7 = 0) 
channel, which has a stronger attractive potential than the 7=1 chan¬ 
nels. We are currently in the process of studying corresponding Schrodinger 
equations for all these channels. 
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